Cluster properties from two-particle angular correlations in p+p collisions at 

= 200 and 410 GeV 
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We present results on two-particle angular correlations in proton-proton collisions at center 
of mass energies of 200 and 410 GeV. The PHOBOS experiment at the Relativistic Heavy Ion 
Collider has a uniquely large coverage for charged particles, giving the opportunity to explore the 
correlations at both short- and long-range scales. At both energies, a complex two-dimensional 
correlation structure in Arj and A<f> is observed. In the context of an independent cluster model 
of short-range correlations, the cluster size and its decay width are extracted from the two-particle 
pseudorapidity correlation function and compared with previous measurements in proton-proton 
and proton-antiproton collisions, as well as PYTHIA and HIJING predictions. 



PACS numbers: 25.75.-q,25.75.Dw,25.75.Gz 



I. INTRODUCTION 

Studies of multiparticle correlations have proven to 
be a powerful tool in exploring the underlying mech- 
anism of particle production in high energy hadronic 
collisions. Inclusive two-particle correlations have been 
found (e.g. Ref. to have two components: direct 
two-particle correlations as well as an effective "long- 
range" correlation due to event-by-event fluctuations of 
the overall particle multiplicity. By considering the 
two-particle rapidity density at fixed multiplicity, the 
intrinsic correlations between particles were isolated and 
found to be approximately Gaussian, with a range of 



arj ~ 1 unit in pseudorapidity (77 



•ln(tan(0/2))). 



Thus, these correlations have been conventionally called 
"short-range" correlations. Their properties have been 
explained by the concept of "cluster" emission. 

The simple idea that hadrons are produced in clus- 
ters, rather than individually, has had great success in 
describing many features of particle production [j], 
In a scenario of independent cluster emission, clusters 
are formed before the final-state hadrons and are inde- 
pendently emitted according to a dynamically generated 



distribution in 77 and <fi. The clusters subsequently decay 
isotropically in their own rest frame into the observed 
final-state hadrons. An independent cluster emission 
model has been widely applied to the study of two- 
particle correlations [1,0], where the observed correlation 
strength and extent in phase space can be parameterized 
in terms of the cluster multiplicity, or "size" (the average 
number of particles in a cluster) and the decay "width" 
(the separation of the particles in pseudorapidity). How- 
ever, it should be noted that independent cluster emission 
is only a phenomenological approach which provides no 
insight about cluster production mechanisms. Further 
modeling will be required to connect these studies to the 
underlying QCD dynamics. 

A related measurement, studying the forward- 
backward multiplicity correlations, was performed in 
Au+Au collisions at center of mass energy per nucleon 
pair (yfisNN) of 200 GeV using the PHOBOS detector 
at the Relativistic Heavy Ion Collider (RHIC ) [|. The 
event-by-event observable C — (Np — Nb) / V Np + Nb 
is constructed, where Np and Nb are defined to be the 
total multiplicity in two symmetric regions forward and 
backward of mid-rapidity. 



The variance (cr^), which is 
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related to the cluster size, was measured. An effective 
cluster size of approximately 2-3, increasing with the 
size of the pseudorapidity window, was observed. In 
heavy ion collisions at RHIC, it has been predicted that 
the formation of a quark gluon plasma (QGP) could 
modify cluster properties relative to P+P collisions Q. 
Unfortunately, the method used in Ref. Q was found to 
have intrinsic limitations for measuring the properties of 
clusters emitted near mid-rapidity. This makes direct 
tests of these predictions difficult and suggests a need for 
different methods to access cluster properties directly. 

In order to explore the properties of short-range 
correlations in more detail and to study cluster prop- 
erties quantitatively, this paper presents the two-particle 
angular correlations in p+p collisions at center of mass 
energies (\/s) of 200 and 410 GeV over a very broad 
acceptance in A77 (= 771 — 772) and A<p (= <pi — </>a). The 
PHOBOS Octagon detector, covering pseudorapidity 
—3 < 77 < 3 over almost the full azimuth, is well suited 
to measure the correlations between particles emitted 
from clusters, which have been found to extend up to 
A77 = 2-3 units of pseudorapidity. The observed two- 
dimensional (2-D) correlation function shows a complex 
correlation structure in the hadronic final state. In 
this analysis, both the cluster size and its decay width 
are extracted from the two-particle correlations as a 
function of A77. Furthermore, the energy, multiplicity 
and Acf> dependence of the cluster size and decay width 
are studied and compared to previous measurements 
as well as PYTHIA @ and HIJING Q predictions at 
various energies. Our present study in p+p collisions will 
provide a useful baseline measurement for understanding 
the hadronization stage in A+A collisions. 



II. DATA SETS 

The data presented here for p+p collisions at y/s = 200 
and 410 GeV were collected during RHIC Run 4 (2004) 
and Run 5 (2005) using the large-acceptance PHOBOS 
Octagon multiplicity detector covering pseudorapidity 
— 3 < 77 < 3 over almost full azimuth. A full description 
of the PHOBOS detector can be found in Ref. 0. Single 
diffractive p+p events were suppressed by requiring at 
least one hit in each of two sets of 16 scintillator trigger 
counters located at distances of —3.21 m and 3.21 m from 
the nominal interaction point z v tx = along the beam 
axis. They cover an acceptance of 3 < \r]\ < 4.5 and 
-180° < < 180°. About 0.5 million 200 GeV and 
0.8 million 410 GeV p+p events were selected for further 
analysis by requiring that the main collision vertex fell 
within |z vtx | < 10 cm along the beam axis. 

The angular coordinates (77, (f>) of charged particles 



are measured using the energy deposited in the silicon 
pads of the Octagon. The granularity of the Octagon is 
determined by the sizes of the readout pads which are 
about 11.25° (~ 0.2 radians) in tj> and range from 0.006 
to 0.05 in 77. Monte-Carlo studies using the PYTHIA 
event generator show that, in the low-multiplicity p+p 
environment, 98% of hits on adjacent pads in 77 are 
created by single primary particles with a small angle 
of incidence. All such neighboring hits are merged to 
single hits in order to avoid fake particle pairs. Noise and 
background hits are rejected by cutting on the deposited 
energy corrected for the path length through the silicon 
after hit merging, assuming that the charged particle 
originated from the main vertex. Depending on 77, the 
merged hits with less than 50-60% of the energy loss 
expected for a minimum ionizing particle are rejected. 
More details of the hit reconstruction procedure can be 
found in Ref. [l(J. 

III. DEFINITION OF THE TWO-PARTICLE 
CORRELATION FUNCTION 

Following an approach similar to that in Ref. Q, the 
inclusive charged two-particle correlation function in two- 
particle (771, 772, 4>\, 4>2) space is defined as follows: 

R(vi,V2, <h.,<h) = {{n-i)[ — — — — — -— - 1 > 

(1) 

where, 

1 / / \ 1 d 2 <r n 

PniV,(P) = -7-7- 

na n a/qdep 

is the single charged particle density distribution and 

Pn(Vl,172,<t>U<h) = —, l —T, , f "77 T7 

n(n — l)<j n dr\\dr}id<p\dq>2 

denotes the charged particle pair distribution. Here a n 
is the cross section of observing n charged particles. The 
above distributions obey the normalization relations: 

J P l n (V,<t>)dr}d<f) = 1 

and 

J Pniviim^it^drjid^d^id^ = 1. 

In practice, we concentrate on the difference in az- 
imuthal angle and pseudorapidity between two particles. 
The correlation function in Eq. [T] is simplified by av- 
eraging over the full acceptance —3 < 771,772 < 3 and 



3 



— 180° < 0i,02 < 180°, reducing the dimensionality 
of the parameter space to Arj and A(f> with a range of 
|At?| < 6 and \A<j>\ < 180°: 

p^iA^A® 



R(Ar),A((>) = ( (n-1 



p mixed (A?7, A0) 



- 1 



(2) 



The pair distribution p„(A?j, A(f>) (with unit integral) 
is determined by taking particle pairs from the same 
event, then averaging over all events. The mixed- 
event background, p mlxed (Ary, A<j)) (with unit integral), 
is constructed by randomly selecting single particles 
from two different events with similar vertex position, 
representing a product of two single particle distribu- 
tions. The vertex bin size of 0.5 cm used in event- 
mixing is mainly determined by the vertex resolution of 
the detector in p+p collisions. Since the background is 
found to be multiplicity independent, we use the inclusive 
p mixcd (A?7, A(j>) in our calculations of Eq. EH The track 
multiplicity n is introduced to compensate for trivial 
dilution effects from uncorrelated particles [2j. This 
stems from the fact that the number of uncorrelated 
pairs grows quadratically with n, while the number of 
correlated pairs grows only linearly. 



IV. CORRECTIONS AND SYSTEM ATICS 

The PHOBOS Octagon detector has a single layer 
of silicon, and so does not provide momentum, charge 
or mass information for the observed particles. Thus 
secondary effects, such as <5-electrons, 7 conversions and 
weak decays cannot be directly rejected. These will 
contribute to correlations unrelated to those between pri- 
mary hadrons, and modify the shape of the measured cor- 
relation function. The incomplete azimuthal acceptance 
in some pseudorapidity regions naturally suppresses the 
overall correlation strength, but MC simulations show 
that it does not change the shape of the correlation in 
Arj and A<f>. To correct for these detector effects in the 
data, correlation functions are calculated for PYTHIA 
events at y/s = 200 GeV both at the generator level for 
true primary charged hadrons, R^f(Ar), A(f>) (Fig. Q}t) 
and with the full GEANT detector simulation and recon- 
struction procedure, R^(A-q, A</>) (Fig. [T]d). As will be 
illustrated in Fig.[3J the actual information in R(Ai], Acf>) 
is reflected to the full A77 and A<j) range in order to more 
clearly show the shape of the correlation function. The 
whole correction procedure can be summarized by the 
following equation: 

iO(A?7,A0) = AxlC^.A^-^A^)]. (3) 

The overall correlation structure consists of both in- 
trinsic and secondary correlations and these two sources 



of correlations are found to be largely independent of 
each other in MC studies, i.e. the correlation from 
secondaries is mostly determined by sensor thickness, 
detector geometry, known cross-sections and decay kine- 
matics. First, within a narrow vertex range (2cm), we 
compare the generator level MC correlation function 
excluding particles outside the PHOBOS detector ac- 
ceptance, i?^.p acc (A?7, Acj>), to the correlation function 
observed after processing the same MC events with all 
the primary hadrons through the GEANT simulation, 
R^(A-q,A(f>) (Fig. [lb). The difference between the two 
correlation functions, S(Arj,A<j)) (Fig. [5b), is attributed 
to the effects of secondary interactions, weak decays, and 
the reconstruction procedure: 

S(A V) A</>) = R^(A V , A<j>) - i$P acc (Ay7, A0). (4) 

To estimate the suppression of the correlation strength 
due to limited acceptance, we compare R^(Ar], A</>) 
and R^ acc (Ar], A</>) using a \ 2 test to extract a scaling 
factor, A, for each vertex bin: 

2 f [iffrFCAg, A0) - A x R^jAru A0)] 2 
X ~ J <P(A,, A(j>) x Ax <P a cc(^, A 0) dA ^ A0 

(5) 

where <7^:p(A?7, A(f>) and CT^p acc (Ary, A</>) are the uncer- 
tainties of the correlation function i?^.p(Ar;, A(f) and 
fl^f acc {Ar], A0). This scaling factor, A, is independent 
of Arj and A4>, and is found to vary between 1.3 and 1.5 
over the vertex range used in this analysis. 

To get the final correlation function, i? d ^(A?7, A<p), 
we apply Eq. [3] to the raw correlation function, 
Rf^ v a (Arj, A(f)). The correction term S{Arj,A4>) from 
Eq. 0] is used to subtract the effects of secondaries and 
weak decays, and the scaling factor, A, from the \ 2 test 
in Eq. [5] is used to remove the suppression due to the 
holes in the PHOBOS acceptance. This procedure is done 
separately for each vertex bin. Three different MC gener- 
ators are used to estimate the systematic uncertainties to 
the correlation function from this correction procedure, 
including PYTHIA, HIJING and a modified PYTHIA 
in which all intrinsic correlations have been removed by 
performing event mixing at the primary hadron level 
(i.e. before weak decays). Typically the systematic error 
(biggest for the peak at Ary=0 and A(f>=0) is less than 
5%. 

Both the correlation function for the data (Fig. [2^) 
and for the PYTHIA simulations (Fig. [TJa) have a sharp 
peak at small Arj and A<p, which is twice as high in the 
data as in the fully simulated MC events, but which 
is not present in the analysis using primary particles 
from the MC generator (Fig. [TJa,) . In the full GEANT 
simulations, the particles contributing to this peak are 
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FIG. 1: Two-particle correlation function in Ar) and A<f> for a) generator level PYTHIA for primary hadrons and b) PYTHIA 
with full GEANT detector simulation and reconstruction procedures in p+p collisions at i/s = 200 GeV. 



a) raw data 200 GeV 




b) PYTHIA 200 GeV 




c) final data 200 GeV 




FIG. 2: a) Raw two-particle correlation function for p+p collision data, b) correction function S(Ar], A(j>) from PYTHIA and 
c) final two-particle correlation function for p+p collision data in A77 and A<f> at y^i = 200 GeV. 



found to be mainly <5-electrons (50%) and 7 conversions 
(40%). The width of the peak is about 0.3 in Att and 28° 
in A(f> for both data and MC, indicating the same origin. 
However, the final data still contain a much narrower 
peak at the near-side of R^(Ar), A<j>) (Fig. [2fc). It 
is likely that this small angle structure results from 
background and detector effects which are not included 
in the MC simulation, although it is not possible to 
rule out unknown physics effects not implemented in 
the event generators. Since the physics of the cluster- 
like particle production investigated in this analysis is 
dominated by correlations on scales of approximately 



one unit in Arj, as will be shown later, we proceed 
by rejecting pairs in a small two-particle acceptance of 
I At; I < 0.15 and \A<f)\ < 5.625° (the single bin centered 
at At/=0 and A<p—0). Studies using primary particles 
from MC generators and the fully simulated events show 
that the extracted cluster parameters, described in the 
next section, change by less than 0.1% due to this cut. 

In addition to the systematic errors related to the 
correction procedure, other systematic uncertainties are 
calculated by varying the vertex position and hit thresh- 
old cuts, and by studying the time dependence of the 
results within the two data sets. Since the acceptance of 





FIG. 3: Two-particle angular correlation function in Aij and A<f> in p+p collisions at y/s — 200 GeV (column a) and 410 GeV 
(column b) (regions of \Ar)\ < 0.15 and \A<j>\ < 5.625° are excluded) shown in two types of representations. The areas enclosed 
by the dashed lines on the lower two panels indicate the actual information in the data which is then reflected to other regions 
to more clearly show the full shape of the correlation. 



our detector is strongly dependent on vertex position, 
any systematic uncertainties due to acceptance and 
geometrical description should manifest themselves as a 
dependence of the final results on vertex position. The hit 
threshold cuts are used to discriminate between primary 
particles and noise or background hits; a variation of 
the thresholds allows us to test the impact of noise and 
background hits on the final results. For each selection of 
reconstruction parameters, a set of correlation functions 
are constructed under different conditions (i.e. 10 differ- 



ent vertex positions), and fully corrected using the MC 
procedure described above. The RMS of these correlation 
functions is calculated as an estimate of the systematic 
error from each particular source. The systematic errors 
due to the vertex position dependency, hit threshold cuts 
and variations in different subsets of events are added 
in quadrature to get the total RMS for each bin in Arj 
and A(f>. The vertex position dependency turns out to 
be the dominant source of the uncertainties. The final 
systematic uncertainties are quoted as 90% C.L.(1.6 x 



6 



RMS). 



V. RESULTS 

The final two-particle inclusive correlation functions, 
averaged over 10 vertex bins, are shown in Fig. [3] as a 
function of Arj and A<f> at y/s = 200 GcV (column a) and 
410 GeV (column b). The near-side hole corresponds to 
the excluded region of |Ary| < 0.15 and |A0| < 5.625°. 
The systematic uncertainties in the absolute value of 
R(Ar?, A<fi) are of the order of 0.3, relative to a peak value 
of 5, with little A77 or A<j> dependence. 

The complex 2-D correlation structure shown in Fig. [3] 
is approximately Gaussian in Arj and persists over the full 
Acj) range, becoming broader toward larger A(f> (which 
will be discussed in quantitative detail below). Similar 
structures also exist in PYTHIA (Fig. [Th) though they 
do not reproduce the strength of the short-range rapidity 
correlations seen in the data. The qualitative features of 
the observed correlation structure are consistent with an 
independent cluster approach according to a simulation 
study from the ISR experiment using a low-mass reso- 
nance (p,uj,r]) gas model [2j. The excess of the near-side 
peak (Ar? ~ and A<j> ~ 0) relative to the away-side 
could be partially a result of the HBT effect [ll|. This 
possibility is investigated in the Appendix using a simple 
MC model and found to be negligible for the cluster 
properties investigated below. 

To study the correlation structure quantitatively, the 
2-D correlation function is projected into a 1-D pseu- 
dorapidity correlation function of Ar] by integrating 
p^(Ar?,A0) and p roixed (Ar?, A</>) over A0 as follows: 



R(Ari) = ( (n - 1) 



J> mixcd (Ar?, Acj))dA(j) 



- 1 



(6) 

The two-particle pseudorapidity correlation function 
R(Ar?), averaged over the Acj) range from 0° to 180°, 
is shown in Fig. @] at yfs = 200 and 410 GeV. The error 
bars (also in Figs. [5] - [5]) correspond to point-to-point 
systematic errors with 90% C.L. The error bands (also in 
Figs. [5]-[5]) denote an overall scale error with 90% C.L. 
as an indication of the uncertainties in the correction 
method which tends to move all of the data points up 
and down in a correlated fashion. The statistical errors 
are negligible due to the large p+p event sample used in 
this analysis. 

In the context of an independent cluster emission 
model, R(Arj) takes the functional form [4j: 



i?(Ar?) = a 



T(Ar?) 



l(Ar?) 



1 



(7) 



n — 1 — 1 — 1 — 1 — r 
p+p 200 GeV 




FIG. 4: Two-particle pseudorapidity correlation function, 
averaged over the A<j> range from 0° to 180° , in p+p collisions 
at = 200 GeV (left) and 410 GeV (right). The solid curves 
correspond to the fits by the cluster model using Eq.[7]over the 
full Arj range. The error bars and bands correspond to point- 
to-point systematic errors and overall scale errors respectively 
with 90% C.L. The statistical errors are negligible. 



(K(K-l)) 
(K) 



is a 



where the correlation strength a 
parameter containing information on the distribution 
of cluster size K. The function r(Ar?) is a Gaussian 
function 



ocexp[-(Ar?)7(4<5 2 )] 

characterizing the correlation of particles originating 
from a single cluster where 5 indicates the decay width 
of the clusters. The background distribution p mlxod (Ar?) 
is just the distribution obtained by event-mixing intro- 
duced in section III. To correct for the holes in the PHO- 
BOS acceptance, we calculate the ratio of the background 
for PYTHIA primary particles, p^p Ti (Ar]), to the one 
obtained in the full GEANT simulations, p^f im (Aj]). 
The ratio is applied to the background calculated from 
the data, p™ata raw( Ar 7)' as a multiplicative factor: 



Pdata, final 

(Ar?) 



-mixed 
A'MCsin 



(Ar?) 



mixed 
rM at a, raw 



(Ar?). (8) 



The effective cluster size is related to the extracted 
correlation strength via the relation: 



K. 



off 



1 



(K(K-l)) 
(K) 



+ 1 = (K} + 



(9) 



Without any knowledge of the distribution of K, it is 
impossible to directly measure the average cluster size 
(K). However, by a x 2 fit of Eq. [7] to the measured two- 
particle pseudorapidity correlation function, the effective 
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FIG. 5: K e g (top) and 8 (bottom) as a function of y/s 
measured by PHOBOS in solid circles, as well as UA5 [l| 
(solid squares) and ISR 0, (solid stars) experiments for 
p+p and p+p collisions. Open circles and squares show 
the PYTHIA and HIJING results respectively. The error 
representations are identical to those in Fig. 2] 

cluster size K e g and decay width 5 can be estimated. 
The independent cluster model provides a good fit to 
the data over a large range in A77, as shown in Fig. 3] 
An effective cluster size K e g = 2.44 + 0.08 and width 
6 = 0.66±0.03 for yfi = 200 GeV and K eS = 2.66 + 0.10, 
6 = 0.67+0.03 for y/s — 410 GeV are obtained with scale 
errors of 5% for K e g and 3% for 6. The three most central 
points always lie above the fits, which could be due to 
residual secondary effects which are left uncorrected or 
other physics at this small scale in Aij. All fit results 
include these three points, but excluding them from the 



fit affects K c g and 5 by no more than 3%. 

In Fig. O our data are compared with previous mea- 
surements of K e g and S as a function of y/s. At lower ISR 
energies (HH)], K e g is constant within error bars. At the 
higher SPS energies [![, UA5 finds K e g to be larger than 
at the ISR, but with little energy dependence between 
200 GeV and 900 GeV. The PHOBOS data are in good 
agreement with the UA5 measurements and, with much 
higher statistics in the p+p event sample, show a clear 
energy dependence of K c g . By contrast, the cluster decay 
width S remains almost constant over the full range of 
collision energies. The event generators, HIJING and 
PYTHIA, show a similar energy dependence of K e g and 
6 to the data, but with a significantly lower magnitude 
of K cS . 

The observed cluster size cannot be fully explained 
by a resonance decay model even at very low energies, 
since the expectation of (K) from resonance decays is 
about 1.5 (extrapolating to 1.7 for K e g depending on the 
assumed cluster size distribution [l[ ) . This is significantly 
lower than the observed values, but is close to what is 
seen in PYTHIA. The HBT effect, after averaging over 
A(f>, would increase the cluster size by no more than 
2% (see Appendix). Additional sources of short-range 
correlations, such as local quantum number conservation 
[HI], are needed to describe the data. As the energy 
increases, the onset of jets should play a more important 
role in the particle production giving bigger clusters, 
which could be the underlying cause for the observed 
energy dependence of K e g. At the LHC, with p+p 
collision at = 14 TeV, jet-like particle production is 
expected to be dominant and should manifest itself in a 
further increase in the effective cluster size. 

To gain further detailed information, the normalized 
cluster parameters, K e g/(K e g) and 5/(5), are calculated 
as a function of the normalized charged multiplicity 
n/{n) at y/s = 200 and 410 GeV (Fig. [BJ>. Scaled by 
the average charged multiplicity (n), the distribution 
of n/(n) turns out to be essentially identical for data 
and charged primary tracks from MC (PYTHIA and 
HIJING), despite the holes in the PHOBOS Octagon 
detector. By dividing the cluster parameters at different 
multiplicity by the averaged values, a scaling behavior is 
observed between two different energies in both data and 
MC. K c g/ (K c g) increases with event multiplicity while 
6/(5) is found to be largely independent of it. PYTHIA 
and HIJING give a similar multiplicity dependence, but 
the increase in K c g/(K e g) is not as strong as in the 
data. Measurements from the ISR and UA5 experiments 
[Ell [HQ are qualitatively consistent with PHOBOS, 
but significantly limited by statistics, and are thus not 
shown here. 

To explicitly show the Acj> dependence of the short- 
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FIG. 6: Normalized effective cluster size K e s/(K e s) (top) 
and decay width 8/(5) (bottom) as a function of normalized 
multiplicity n/ (n) in p+p collisions at y/s = 200 and 410 GeV 
measured by PHOBOS (solid symbols), as well as MC studies 
(open symbols). The error representations are identical to 
those in Fig. 0] 

range pseudorapidity correlation seen in Fig. [3J the 
A(p range from 0° to 180° is divided into four regions 
and projected separately onto the Arj axis. Fig. [7] 
shows K c s/(K c fi) and 8/(8) for different A<fi regions. 
K e s/(K c s) gradually decreases while 5 increases as 
one goes from small to large A</> region. This might 
reflect some information about the transverse momentum 
distribution of the clusters 01 . High pr clusters should 
generally contribute to a narrow hump in the near-side 
(near Acf> = 0°) of the correlation function in Fig. [3l 
whereas the broader away-side (near A<f) = 180°) comes 



FIG. 7: Normalized effective cluster size K e ff/(K e f{) (top) 
and decay width 5/(5) (bottom) as a function of A<f> in p+p 
collisions at -Js = 200 and 410 GeV measured by PHOBOS 
(solid symbols), as well as MC studies (open symbols). The 
error representations are identical to those in Fig. [4] 

from clusters with lower transverse momentum. Again, 
PYTHIA and HIJING are qualitatively similar to data 
but show a smaller decay width at intermediate A(j) and 
a larger decay width near A(f> = 180°. 

The two-particle azimuthal correlation functions 
R(A0), averaged over a broad range of A77 from to 6, 
in p+p collisions at y/s = 200 and 410 GeV are presented 
in Fig.[H] R(A(j)) is obtained using a procedure similar to 
that shown for R(A^) in Eq.[6] The observed asymmetric 
structure in Acj) (with the Ar^-averaged away-side peak 
larger than the near-side peak) could also provide some 
information about the momentum and size distribution 
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collisions at RHIC energies. 

We acknowledge the generous support of the Collider- 
Accelerator Department. This work was partially 
supported by U.S. DOE grants DE-AC02-98CH10886, 
DE-FG02-93ER40802, DE-FC02-94ER40818, DE-FG02- 
94ER40865, DE-FG02-99ER41099, and W-31-109-ENG- 
38, by U.S. NSF grants 9603486, 0072204, and 0245011, 
by Polish KBN grant l-P03B-062-27(2004-2007), by NSC 
of Taiwan Contract NSC 89-2112-M-008-024, and by 
Hungarian OTKA grant (F 049823). 



Appendix: HBT Correlations 



FIG. 8: Two-particle azimuthal correlation function, aver- 
aged over the Ar/ range from to 6, in p+p collisions at 
= 200 GeV (left) and 410 GeV (right) from PHOBOS 
(solid circles), as well as MC studies (open circles). The error 
representations are identical to those in Fig. [4] 

of clusters. The A77 integrated correlation function is 
similar in magnitude both for data and PYTHIA, despite 
the significant difference in the extracted K c g. More 
detailed modeling of the cluster properties is needed to 
fully explain many aspects of the complex two-particle 
correlation function. 



VI. CONCLUSION 

PHOBOS has measured two-particle angular correla- 
tions over a wide range in Arj and A</> in p+p collisions 
at y/s = 200 and 410 GeV. Short-range correlations are 
observed over the full range in A<f>, with a maximum 
at A77 = which becomes wider at larger Acf>. In 
the context of the cluster model, the effective cluster 
size and decay width are extracted from the two-particle 
pseudorapidity correlation function, and compared with 
previous experiments, as well as HIJING and PYTHIA 
event generators. Dependence of the cluster size on both 
beam energy and scaled multiplicity is observed, while 
the cluster width is essentially constant. The short-range 
correlation strength (or equivalently the effective cluster 
size, K e g) exceeds the expectation from the decays 
of resonance particles, suggesting the need for other 
sources of short-range correlations. As mentioned in 
the introduction, cluster properties could be modified in 
A+A collisions relative to p+p collisions by the presence 
of a QGP[(|. Future studies should clarify the evolution 
of cluster parameters from p+p to d+Au and A+A 



The Bose-Einstein or HBT correlation will certainly 
contribute to the two-particle angular correlation func- 
tion and thus might impact the estimation of cluster 
parameters (11 ] . It is known to cause a strong correlation 
between two identical particles having small relative in- 
variant four-momentum q$ nv (=(qi~ q2) 2 )- For identical 
bosons, the overall shape of the HBT effect can be 
approximated by a 1-D Gaussian correlation function: 

C(q in J = l + Ae-^«?-. (10) 

To estimate the effect of the HBT correlation, in 
PYTHIA, which does not include this effect, we calculate 
the q in „ for each pair in f^(Ar], Acf)) (ignoring the 
particle species), and weight the pair by C(q inv ) to 
artificially introduce the HBT correlation. In our study, 
we use A = 0.8 and Ri nv = 1.0 fm. Considering the 
fact that most of the pairs measured in this analysis 
are not identical-particle pairs, which should reduce the 
magnitude of A, this choice of parameters exaggerates the 
expected effect of the HBT correlation. Thus, it gives a 
conservative estimate of its potential effect on the data. 

In Fig. [HI the difference between the correlation func- 
tion with and without C(q jnlJ ) weighting is shown as a 
function of A?y and A<j> for PYTHIA. Finally, Fig. [TO] 
shows a comparison of pseudorapidity correlation func- 
tions with and without HBT weighting, with extracted 
cluster sizes of 1.76 and 1.72 respectively. From these 
studies, it appears that the HBT effect does enhance the 
short-range correlations with a range of around one unit 
in Ar] and 45° in A<p. However, after averaging over 
Acjj, it only slightly influences the rapidity correlation 
function, resulting in an increase in cluster size of at most 
2%. Even in the region 0° < A<f> < 45°, the increase is 
less than 3.5%. 
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FIG. 9: Difference of the correlation function with HBT 
weighting minus the one without HBT weighting in A?7 and 
A<p for PYTHIA at y/s = 200 GeV 



; • Standard PYTHIA 

■_ o PYTHIA with HBT weighting 




FIG. 10: Comparison of the pseudorapidity correlation func- 
tion with (open circles) and without (solid circles) HBT 
weighting for PYTHIA at y/s = 200 GeV. 
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